COVID-19 outbreak at a residential apartment building in Northern Ontario, Canada.

In February 2021, a cluster of Beta variant (B.1.351) coronavirus disease 2019 (COVID-19) cases were identified in an apartment building located in Northern Ontario, Canada. Most cases had no known contact with each other. Objectives of this multi-component outbreak investigation were to better understand the social and environmental factors that facilitated the transmission of COVID-19 through this multi-unit residential building (MURB). A case-control study examined building-specific exposures and resident behaviours that may have increased the odds of being a case. A professional engineer assessed the building's heating, ventilation, and air-conditioning (HVAC) systems. Whole-genome sequencing and an in-depth genomic analysis were performed. Forty-five outbreak-confirmed cases were identified. From the case-control study, being on the upper floors (OR: 10.4; 95% CI: 1.63-66.9) and within three adjacent vertical lines (OR: 28.3; 3.57-225) were both significantly associated with being a case of COVID-19, after adjusting for age. There were no significant differences in reported behaviours, use of shared spaces, or precautions taken between cases and controls. An assessment of the building's ventilation found uncontrolled air leakage between apartment units. A single genomic cluster was identified, where most sequences were identical to one another. Findings from the multiple components of this investigation are suggestive of aerosol transmission between units.

Putting everything in its place: using the INSDC compliant Pathogen Data Object Model to better structure genomic data submitted for public health applications.

Fast, efficient public health actions require well-organized and coordinated systems that can supply timely and accurate knowledge. Public databases of pathogen genomic data, such as the International Nucleotide Sequence Database Collaboration (INSDC), have become essential tools for efficient public health decisions. However, these international resources began primarily for academic purposes, rather than for surveillance or interventions. Now, queries need to access not only the whole genomes of multiple pathogens but also make connections using robust contextual metadata to identify issues of public health relevance. Databases that over time developed a patchwork of submission formats and requirements need to be consistently organized and coordinated internationally to allow effective searches.To help resolve these issues, we propose a common pathogen data structure called the Pathogen Data Object Model (DOM) that will formalize the minimum pieces of sequence data and contextual data necessary for general public health uses, while recognizing that submitters will likely withhold a wide range of non-public contextual data. Further, we propose contributors use the Pathogen DOM for all pathogen submissions (bacterial, viral, fungal, and parasites), which will simplify data submissions and provide a consistent and transparent data structure for downstream data analyses. We also highlight how improved submission tools can support the Pathogen DOM, offering users additional easy-to-use methods to ensure this structure is followed.

Repeated Emergence of Variant TetR Family Regulator, FarR, and Increased Resistance to Antimicrobial Unsaturated Fatty Acid among Clonal Complex 5 Methicillin-Resistant Staphylococcus aureus.

Resistance-nodulation-division (RND) superfamily efflux pumps promote antibiotic resistance in Gram-negative pathogens, but their role in Gram-positive pathogens, including methicillin-resistant Staphylococcus aureus (MRSA) is undocumented. However, recent in vitro selections for resistance of S. aureus to an antimicrobial fatty acid, linoleic acid, and an antibiotic, rhodomyrtone, identified H121Y and C116R substitution variants, respectively, in a TetR family regulator, FarR, promoting increased expression of the RND pump FarE. Hypothesizing that in vivo selection pressures have also promoted the emergence of FarR variants, we searched available genome data and found that strains with FarRH121Y from human and bovine hosts have emerged sporadically in clonal complexes (CCs) CC1, CC30, CC8, CC22, and CC97, whereas multiple FarR variants have occurred within CC5 hospital-associated (HA)-MRSA. Of these, FarRE160G and FarRE93EE were exclusive to CC5, while FarRC116Y, FarRP165L, and FarRG166D also occurred in nonrelated CCs, primarily from bovine hosts. Within CC5, FarRC116Y and FarRG166D strains were polyphyletic, each exhibiting two emergence events. FarRC116Y and FarRE160G were individually sufficient to confer increased expression of FarE and enhanced resistance to linoleic acid (LA). Isolates with FarRE93EE were most closely related to S. aureus N315 MRSA and exhibited increased resistance independently of FarRE93EE. Accumulation of pseudogenes and additional polymorphisms in FarRE93EE strains contributed to a multiresistance phenotype which included fosfomycin and fusidic acid resistance in addition to increased linoleic acid resistance. These findings underscore the remarkable adaptive capacity of CC5 MRSA, which includes the polyphyletic USA100 lineage of HA-MRSA that is endemic in the Western hemisphere and known for the acquisition of multiple resistance phenotypes.

The DataHarmonizer: a tool for faster data harmonization, validation, aggregation and analysis of pathogen genomics contextual information.

Pathogen genomics is a critical tool for public health surveillance, infection control, outbreak investigations as well as research. In order to make use of pathogen genomics data, they must be interpreted using contextual data (metadata). Contextual data include sample metadata, laboratory methods, patient demographics, clinical outcomes and epidemiological information. However, the variability in how contextual information is captured by different authorities and how it is encoded in different databases poses challenges for data interpretation, integration and their use/re-use. The DataHarmonizer is a template-driven spreadsheet application for harmonizing, validating and transforming genomics contextual data into submission-ready formats for public or private repositories. The tool's web browser-based JavaScript environment enables validation and its offline functionality and local installation increases data security. The DataHarmonizer was developed to address the data sharing needs that arose during the COVID-19 pandemic, and was used by members of the Canadian COVID Genomics Network (CanCOGeN) to harmonize SARS-CoV-2 contextual data for national surveillance and for public repository submission. In order to support coordination of international surveillance efforts, we have partnered with the Public Health Alliance for Genomic Epidemiology to also provide a template conforming to its SARS-CoV-2 contextual data specification for use worldwide. Templates are also being developed for One Health and foodborne pathogens. Overall, the DataHarmonizer tool improves the effectiveness and fidelity of contextual data capture as well as its subsequent usability. Harmonization of contextual information across authorities, platforms and systems globally improves interoperability and reusability of data for concerted public health and research initiatives to fight the current pandemic and future public health emergencies. While initially developed for the COVID-19 pandemic, its expansion to other data management applications and pathogens is already underway.

Divergent SARS-CoV-2 variant emerges in white-tailed deer with deer-to-human transmission.

Wildlife reservoirs of broad-host-range viruses have the potential to enable evolution of viral variants that can emerge to infect humans. In North America, there is phylogenomic evidence of continual transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from humans to white-tailed deer (Odocoileus virginianus) through unknown means, but no evidence of transmission from deer to humans. We carried out an observational surveillance study in Ontario, Canada during November and December 2021 (n = 300 deer) and identified a highly divergent lineage of SARS-CoV-2 in white-tailed deer (B.1.641). This lineage is one of the most divergent SARS-CoV-2 lineages identified so far, with 76 mutations (including 37 previously associated with non-human mammalian hosts). From a set of five complete and two partial deer-derived viral genomes we applied phylogenomic, recombination, selection and mutation spectrum analyses, which provided evidence for evolution and transmission in deer and a shared ancestry with mink-derived virus. Our analysis also revealed an epidemiologically linked human infection. Taken together, our findings provide evidence for sustained evolution of SARS-CoV-2 in white-tailed deer and of deer-to-human transmission.

Real-Time RT-PCR Allelic Discrimination Assay for Detection of N501Y Mutation in the Spike Protein of SARS-CoV-2 Associated with B.1.1.7 Variant of Concern.

The N501Y amino acid mutation caused by a single point substitution A23063T in the spike gene of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is possessed by three variants of concern (VOCs), B.1.1.7, B.1.351, and P.1. A rapid screening tool using this mutation is important for surveillance during the coronavirus disease 2019 (COVID-19) pandemic. We developed and validated a single nucleotide polymorphism real-time reverse transcription PCR assay using allelic discrimination of the spike gene N501Y mutation to screen for potential variants of concern and differentiate them from SARS-CoV-2 lineages without the N501Y mutation. A total of 160 clinical specimens positive for SARS-CoV-2 were characterized as mutant (N501Y) or N501 wild type by Sanger sequencing and were subsequently tested with the N501Y single nucleotide polymorphism real-time reverse transcriptase PCR assay. Our assay, compared to Sanger sequencing for single nucleotide polymorphism detection, demonstrated positive percent agreement of 100% for all 57 specimens displaying the N501Y mutation, which were confirmed by Sanger sequencing to be typed as A23063T, including one specimen with mixed signal for wild type and mutant. Negative percent agreement was 100% in all 103 specimens typed as N501 wild type, with A23063 identified as wild type by Sanger sequencing. The identification of circulating SARS-CoV-2 lineages carrying an N501Y mutation is critical for surveillance purposes. Current identification methods rely primarily on Sanger sequencing or whole-genome sequencing, which are time consuming, labor intensive, and costly. The assay described herein is an efficient tool for high-volume specimen screening for SARS-CoV-2 VOCs and for selecting specimens for confirmatory Sanger or whole-genome sequencing. IMPORTANCE During the coronavirus disease 2019 (COVID-19) pandemic, several variants of concern (VOCs) have been detected, for example, B.1.1.7, B.1.351, P.1, and B.1.617.2. The VOCs pose a threat to public health efforts to control the spread of the virus. As such, surveillance and monitoring of these VOCs is of the utmost importance. Our real-time RT-PCR assay helps with surveillance by providing an easy method to quickly survey SARS-CoV-2 specimens for VOCs carrying the N501Y single nucleotide polymorphism (SNP). Samples that test positive for the N501Y mutation in the spike gene with our assay can be sequenced to identify the lineage. Thus, our assay helps to focus surveillance efforts and decrease turnaround times.

Coronavirus Disease 2019 (COVID-19) Outbreak Associated With Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) P.1 Lineage in a Long-Term Care Home After Implementation of a Vaccination Program-Ontario, Canada, April-May 2021.

In a P.1 coronavirus disease 2019 (COVID-19) outbreak in a long-term care home, vaccine effectiveness against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection was 52.5% (95% confidence interval: 26.9%-69.1%) in residents and 66.2% (2.3%-88.3%) in staff. Vaccine effectiveness against severe illness was 78.6% (47.9%-91.2%) in residents. Two of 19 vaccinated resident case patients died. Outbreak management required both vaccination and infection control measures.

Characteristics of SARS-CoV-2 testing for rapid diagnosis of COVID-19 during the initial stages of a global pandemic.

Accurate SARS-CoV-2 diagnosis is essential to guide prevention and control of COVID-19. Here we examine SARS-CoV-2 molecular-based test performance characteristics and summarize case-level data related to COVID-19 diagnosis. From January 11 through April 22, 2020, Public Health Ontario conducted SARS-CoV-2 testing of 86,942 specimens collected from 80,354 individuals, primarily using real-time reverse-transcription polymerase chain reaction (rRT-PCR) methods. We analyzed test results across specimen types and for individuals with multiple same-day and multi-day collected specimens. Nasopharyngeal compared to throat swabs had a higher positivity (8.8% vs. 4.8%) and an adjusted estimate 2.9 Ct lower (SE = 0.5, p<0.001). Same-day specimens showed high concordance (98.8%), and the median Ct of multi-day specimens increased over time. Symptomatic cases had rRT-PCR results with an adjusted estimate 3.0 Ct (SE = 0.5, p<0.001) lower than asymptomatic/pre-symptomatic cases. Overall test sensitivity was 84.6%, with a negative predictive value of 95.5%. Molecular testing is the mainstay of SARS-CoV-2 diagnosis and testing protocols will continue to be dynamic and iteratively modified as more is learned about this emerging pathogen.

Enquête sur une éclosion importante de SRAS-CoV-2 dans un établissement de soins de longue durée au début de la pandémie.

CONTEXTE: Le déploiement de mesures de gestion des éclosions de SRAS-CoV-2 dans les établissements de soins de longue durée en Ontario a permis d'en réduire la fréquence et la gravité. Nous décrivons ici les données épidémiologiques et de laboratoire d'une de ces premières éclosions en Ontario afin de déterminer les facteurs associés à son importance et les impacts des interventions progressives de lutte contre les infections appliquées pendant la durée de l'éclosion. MÉTHODES: Nous avons obtenu du bureau de santé la liste des cas et les données de l'éclosion afin de décrire les cas chez les résidents et le personnel, leur gravité et leur distribution dans le temps et à l'intérieur de l'établissement touché. Quand elles étaient disponibles, nous avons obtenu des données concernant les échantillons soumis au laboratoire de Santé publique Ontario et effectué un séquençage complet et une analyse phylogénétique des échantillons viraux de l'éclosion. RÉSULTATS: Sur les 65 résidents de l'établissement de soins de longue durée, 61 (94 %) ont contracté le SRAS-CoV-2, le taux de létalité étant de 45 % (28/61). Parmi les 67 employés initiaux, 34 (51 %) ont contracté le virus, et aucun n'est décédé. Lorsque l'éclosion a été déclarée, 12 employés, 2 visiteurs et 9 résidents présentaient des symptômes. Parmi les résidents, les cas se trouvaient dans 3 des 4 secteurs de l'établissement. L'analyse phylogénétique a montré une forte similitude des séquences; une seule autre souche de SRAS-CoV-2 génétiquement distincte a été identifiée chez un employé à la troisième semaine de l'éclosion. Après le déploiement de toutes les mesures de gestion de l'éclosion, aucun cas n'a été identifié parmi les 26 nouveaux employés appelés en renfort. INTERPRÉTATION: La propagation rapide et non détectée du virus dans un établissement de soins de longue durée a donné lieu à des taux élevés d'infection chez les résidents et le personnel. L'application progressive de mesures de gestion après le pic de l'éclosion a permis d'éviter la contamination du personnel appelé en renfort et fait désormais partie des politiques à long terme de prévention des éclosions en Ontario.

Temporal Dynamics and Evolution of SARS-CoV-2 Demonstrate the Necessity of Ongoing Viral Genome Sequencing in Ontario, Canada.

Genome-wide variation in SARS-CoV-2 reveals evolution and transmission dynamics which are critical considerations for disease control and prevention decisions. Here, we review estimates of the genome-wide viral mutation rates, summarize current COVID-19 case load in the province of Ontario, Canada (5 January 2021), and analyze published SARS-CoV-2 genomes from Ontario (collected prior to 24 November 2020) to test for more infectious genetic variants or lineages. The reported mutation rate (∼10-6 nucleotide [nt]-1 cycle-1) for SARS-CoV-2 is typical for coronaviruses. Analysis of published SARS-CoV-2 genomes revealed that the G614 spike protein mutation has dominated infections in Ontario and that SARS-CoV-2 lineages present in Ontario have not differed significantly in their rate of spread. These results suggest that the SARS-CoV-2 population circulating in Ontario has not changed significantly to date. However, ongoing genome monitoring is essential for identification of new variants and lineages that may contribute to increased viral transmission.

Investigation of a severe SARS-CoV-2 outbreak in a long-term care home early in the pandemic.

BACKGROUND: The implementation of outbreak management measures has decreased the frequency and severity of SARS-CoV-2 outbreaks in Ontario long-term care homes. We describe the epidemiological and laboratory data from one of the first such outbreaks in Ontario to assess factors associated with its severity, and the impact of progressive interventions for infection control over the course of the outbreak. METHODS: We obtained line list and outbreak data from the public health unit to describe resident and staff cases, severity and distribution of cases over time and within the outbreak facility. Where available, we obtained data on laboratory specimens from the Public Health Ontario Laboratory and performed whole genome sequencing and phylogenetic analysis of viral specimens from the outbreak. RESULTS: Among 65 residents of the long-term care home, 61 (94%) contracted SARS-CoV-2, with a case fatality rate of 45% (28/61). Among 67 initial staff, 34 (51%) contracted the virus and none died. When the outbreak was declared, 12 staff, 2 visitors and 9 residents had symptoms. Resident cases were located in 3 of 4 areas of the home. Phylogenetic analysis showed tight clustering of cases, with only 1 additional strain of genetically distinct SARS-CoV-2 identified from a staff case in the third week of the outbreak. No cases were identified among 26 new staff brought into the home after full outbreak measures were implemented. INTERPRETATION: Rapid and undetected viral spread in a long-term care home led to high rates of infection among residents and staff. Progressive implementation of outbreak measures after the peak of cases prevented subsequent staff cases and are now part of long-term care outbreak policy in Ontario.

Prevalence of Co-Infections with Respiratory Viruses in Individuals Investigated for SARS-CoV-2 in Ontario, Canada.

BACKGROUND: Co-infections of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with respiratory viruses, bacteria and fungi have been reported to cause a wide range of illness. OBJECTIVES: We assess the prevalence of co-infection of SARS-CoV-2 with seasonal respiratory viruses, document the respiratory viruses detected among individuals tested for SARS-CoV-2, and describe characteristics of individuals with respiratory virus co-infection detected. METHODS: Specimens included in this study were submitted as part of routine clinical testing to Public Health Ontario Laboratory from individuals requiring testing for SARS-CoV-2 and/or seasonal respiratory viruses. RESULTS: Co-infection was detected in a smaller proportion (2.5%) of individuals with laboratory confirmed SARS-CoV-2 than those with seasonal respiratory viruses (4.3%); this difference was not significant. Individuals with any respiratory virus co-infection were more likely to be younger than 65 years of age and male than those with single infection. Those with SARS-CoV-2 co-infection manifested mostly mild respiratory symptoms. CONCLUSIONS: Findings of this study may not support routine testing for seasonal respiratory viruses among all individuals tested for SARS-CoV-2, as they were rare during the study period nor associated with severe disease. However, testing for seasonal respiratory viruses should be performed in severely ill individuals, in which detection of other viruses may assist with patient management.

Detection of the novel SARS-CoV-2 European lineage B.1.177 in Ontario, Canada.

Background: Travel-related dissemination of SARS-CoV-2 continues to contribute to the global pandemic. A novel SARS-CoV-2 lineage (B.1.177) reportedly arose in Spain in the summer of 2020, with subsequent spread across Europe linked to travel by infected individuals. Surveillance and monitoring through the use of whole genome sequencing (WGS) offers insights into the global and local movement of pathogens such as SARS-CoV-2 and can detect introductions of novel variants. Methods: We analysed the genomes of SARS-CoV-2 sequenced for surveillance purposes from specimens received by Public Health Ontario (Sept 6 - Oct 10, 2020), collected from individuals in eastern Ontario, which comprised the study sample. Taxonomic lineages were identified using pangolin (v2.08) and phylogenetic analysis incorporated publicly available genomes covering the same time period as the study sample. Epidemiological data collected from laboratory requisitions and standard reportable disease case investigation was integrated into the analysis. Results: Genomic surveillance identified a COVID-19 case with SARS-CoV-2 lineage B.1.177 from an individual in eastern Ontario in late September, 2020. The individual had recently returned from Europe. Genomic analysis with publicly available data indicate the most closely related genomes to this specimen were from Southern Europe. Genomic surveillance did not identify further cases with this lineage. Conclusions: Genomic surveillance allowed for early detection of a novel SARS-CoV-2 lineage in Ontario which was deemed to be travel related. This type of genomic-based surveillance is a key tool to measure the effectiveness of public health measures such as mandatory self-isolation for returned travellers, aimed at preventing onward transmission of newly introduced lineages of SARS-CoV-2.

Using Whole-genome Sequencing to Determine the Timing of Secondary Tuberculosis in British Columbia, Canada.

Combined with epidemiological data, whole-genome sequencing (WGS) can help better resolve individual tuberculosis (TB) transmission events to a degree not possible with traditional genotyping. We combine WGS data with patient-level data to calculate the timing of secondary TB among contacts of people diagnosed with active TB in British Columbia, Canada.

Genomic Epidemiology of Invasive Methicillin-Resistant Staphylococcus aureus Infections Among Hospitalized Individuals in Ontario, Canada.

BACKGROUND: Prevention and control of methicillin-resistant Staphylococcus aureus (MRSA) infections remain challenging. In-depth surveillance integrating patient and isolate data can provide evidence to better inform infection control and public health practice. METHODS: We analyzed MRSA cases diagnosed in 2010 (n = 212) and 2016 (n = 214) by hospitals in Ontario, Canada. Case-level clinical and demographic data were integrated with isolate characteristics, including antimicrobial resistance (AMR), classic genotyping, and whole-genome sequencing results. RESULTS: Community-associated MRSA (epidemiologically defined) increased significantly from 23.6% in 2010 to 43.0% in 2016 (P < .001). The MRSA population structure changed over time, with a 1.5× increase in clonal complex (CC)8 strains and a concomitant decrease in CC5. The clonal shift was reflected in AMR patterns, with a decrease in erythromycin (86.7% to 78.4%, P = .036) and clindamycin resistance (84.3% to 47.9%, P < .001) and a >2-fold increase in fusidic acid resistance (9.0% to 22.5%, P < .001). Isolates within both CC5 and CC8 were relatively genetically diverse. We identified 6 small genomic clusters-3 potentially related to transmission in healthcare settings. CONCLUSIONS: Community-associated MRSA is increasing among hospitalized individuals in Ontario. Clonal shifting from CC5 to CC8 has impacted AMR. We identified a relatively high genetic diversity and limited genomic clustering within these dominant CCs.

Whole genome sequencing of Mycobacterium tuberculosis: current standards and open issues.

Whole genome sequencing (WGS) of Mycobacterium tuberculosis has rapidly progressed from a research tool to a clinical application for the diagnosis and management of tuberculosis and in public health surveillance. This development has been facilitated by drastic drops in cost, advances in technology and concerted efforts to translate sequencing data into actionable information. There is, however, a risk that, in the absence of a consensus and international standards, the widespread use of WGS technology may result in data and processes that lack harmonization, comparability and validation. In this Review, we outline the current landscape of WGS pipelines and applications, and set out best practices for M. tuberculosis WGS, including standards for bioinformatics pipelines, curated repositories of resistance-causing variants, phylogenetic analyses, quality control and standardized reporting.

The problem with defining foreign birth as a risk factor in tuberculosis epidemiology studies.

OBJECTIVE: To examine how stratifying persons born outside Canada according to tuberculosis (TB) incidence in their birth country and other demographic factors refines our understanding of TB epidemiology and local TB transmission. BACKGROUND: Population-level TB surveillance programs and research studies in low incidence settings often report all persons born outside the country in which the study is conducted as "foreign-born"-a single label for a highly diverse population with variable TB risks. This may mask important TB epidemiologic trends and not accurately reflect local transmission patterns. METHODS: We used population-level data from two large cohorts in British Columbia (BC), Canada: an immigration cohort (n = 337,492 permanent residents to BC) and a genotyping cohort (n = 2290 culture-confirmed active TB cases). We stratified active TB case counts, incidence rates, and genotypic clustering (an indicator of TB transmission) in BC by birth country TB incidence, age at immigration, and years since arrival. RESULTS: Persons from high-incidence countries had a 12-fold higher TB incidence than those emigrating from low-incidence settings. Estimates of local transmission, as captured by genotyping, versus reactivation of latent TB infection acquired outside Canada varied when data were stratified by birthplace TB incidence, as did patient-level characteristics of individuals in each group, such as age and years between immigration and diagnosis. CONCLUSION: Categorizing persons beyond simply "foreign-born", particularly in the context of TB epidemiologic and molecular data, is needed for a more accurate understanding of TB rates and patterns of transmission.

Universal genotyping reveals province-level differences in the molecular epidemiology of tuberculosis.

OBJECTIVES: Compare the molecular epidemiology of tuberculosis (TB) between two large Canadian provinces-Ontario and British Columbia (BC)-to identify genotypic clusters within and across both provinces, allowing for an improved understanding of genotype data and providing context to more accurately identify clusters representing local transmission. DESIGN: We compared 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) genotyping for 3,314 Ontario and 1,602 BC clinical Mycobacterium tuberculosis isolates collected from 2008 through 2014. Laboratory data for each isolate was linked to case-level records to obtain clinical and demographic data. RESULTS: The demographic characteristics of persons with TB varied between provinces, most notably in the proportion of persons born outside Canada, which was reflected in the large number of unique genotypes (n = 3,461). The proportion of clustered isolates was significantly higher in BC. Substantial clustering amongst non-Lineage 4 TB strains was observed within and across the provinces. Only two large clusters (≥10 cases/cluster) representing within province transmission had interprovincial genotype matches. CONCLUSION: We recommend expanding analysis of shared genotypes to include neighbouring jurisdictions, and implementing whole genome sequencing to improve identification of TB transmission, recognize outbreaks, and monitor changing trends in TB epidemiology.